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Abstract. We examine the effects of the CuO chains on the density of states in the 
vortex phase in YBa2Cu307, via a calculation based on the tight-binding proximity 
model. In this model, chain superconductivity results from single-electron hopping 
between the intrinsically- normal chains and intrinsically-superconducting Cu02 planes. 
The calculations are based on self-consistent solutions of the Bogolyubov-de Gennes 
equations for a bilayer consisting of a single Cu02 layer and a single CuO chain layer. 
We find that, in addition to the dispersing resonances found in single-layer models, the 
chains introduce a second set of dispersing resonances associated with the induced gap 
in the chain layer. These new resonances are highly anisotropic and distort the vortex 
core shape. 
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1. Introduction 

Because of the availability of large high-quality single crystals, YE^CusOy-^ is one of 
the most widely studied of the high temperature superconductors (HTS). As with other 
HTS, the essential structural components are the two-dimensional (2D) Cu02 layers. In 
these layers, the Cu and O orbitals are partially filled and are therefore conducting. The 
interlayer coupling is generally weak, meaning that the electronic bands contributing at 
the Fermi energy derive most of their weight from individual CuC>2 layers. For this 
reason, the standard model for a generic cuprate HTS consists of a single 2D CuC>2 
layer in isolation. 

YBa 2 Cu307_,5 differs from other cuprate HTS in one significant respect: in addition 
to the CuC>2 layers, there are layers of one- dimensional (ID) chains. Band structure 
calculations [lj suggest that the chains are roughly one-quarter filled, and contribute a 
quasi- ID Fermi surface. Despite this, the chains are nearly always ignored in theories 
of YBa2Cu307_,5, presumably because it has many physical properties that are common 
to all HTS, suggesting a minor role for the chains. In addition, the metallicity of the 
chains has been questioned, either because electrons in ID are localized by arbitrarily 
weak disorder, or because they are unstable towards the formation of insulating charge 
density wave states. 

Early evidence for the metallicity of the chains came from measurements of the 
resistivity anisotropy [2], which found p a / Pb ~ 2 at 100 K, where p a and pi, are the 
resistivities parallel to the two axes of the Cu02 planes, and where b is also parallel to 
the chains. Such an anisotropy is expected if the chains are conducting and carry current 
in parallel with the planes. Since O-vacancies, which are prevalent in YBa2Cu307_,5, 
would strongly localize ID electrons, the observed anisotropy appears to indicate a 
relatively large hybridization of chain and plane states. Further support for this comes 
from c-axis (perpendicular to the planes) resistivity measurements [3], H] which find an 
anisotropy p c / p a ~ 50 at 100 K in optimally-doped YBa 2 Cu30 6 .93. In a tight binding 
model, this corresponds to t±/t ~ 7 where t and t±_ are the in-plane and perpendicular 
hopping amplitudes. This anisotropy is smaller than in other HTS and suggests coherent 
c-axis transport [5] . 

There has been some debate as to whether the CuO chains undergo a Peierls 
transition to a charge density wave state. Scanning tunneling microscopy (STM) 
experiments on chain-terminated surfaces found charge modulations which were 
interpreted as charge density waves. [6] However, later experiments [7] found that the 
modulation wavelength depends on the bias voltage of the STM tip, which is a 
characteristic of Friedel oscillations (i.e. standing wave patters produced by impurity 
scattering of itinerant electrons) rather than charge density waves. These later 
experiments are consistent with metallic chains. 

Direct evidence for a chain Fermi surface has recently been found in a number 
of angle resolved photoemission experiments (HJ El HQ]. The Fermi surface appears 
to be consistent with that predicted by band structure calculations, although a 
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complete characterization of the band structure is rendered difficult by the existence 
of surface states, which dominate surface-sensitive experiments such as tunneling and 
photoemission. The correspondence between the surface states and electronic states in 
the bulk is not yet established. 

Complementary information on the chain states comes from penetration depth 
anisotropy measurements [TT1 [T2] which indicate that the chains are superconducting 
below the bulk critical temperature T c , and have a substantial superffuid density. 
Because the chains and Cu02 planes are structurally different, the apparent similarity 
of their superconducting states, as measured by the penetration depth, was puzzling 
for many years. The chains, being far from half-filled, should not have an intrinsic 
pairing interaction capable of producing high critical temperatures. A more likely 
scenario is that the chains derive their superconductivity from proximity coupling to the 
Cu02 planes. [13] Calculations showed that proximity models for chain superconductivity 
introduce small energy scales related to the induced gap in the chains. To date, these 
small scales have not been observed in penetration depth measurements. [14J It was 
later shown[15j that the smallest energy scale, which is relevant to low temperature 
measurements, comes from a subset of chain-derived electronic states that are weakly 
hybridized with the Cu02 planes. Since these particular states have a predominantly 
ID character, they are are strongly affected by localization corrections due to chain 
disorder, and are therefore difficult to detect via the electrodynamic response. 

Small energy scales have been seen in other experiments. For example, it has been 
foundpiH [T7] that the vortex cores in YE^CusOy-^, as measured by /iSR experiments, 
are much larger at small magnetic fields B than expected from the measured upper 
critical field, and that the cores contract rapidly with increasing B below a crossover field 
B* . This anomalous behaviour has been explained by the presence of CuO chains[18, 19], 
and it can be shown that B* ~ E^. It has also been suggested that a similar small energy 
scale seen in tunneling experiments |20l |2T| 122] can be attributed to the chains. 

In summary, we argue that there is reasonable evidence that the CuO chains in 
YBa 2 Cu307_5 (a) are metallic, (b) are signifcantly hybridized with the Cu02 planes, 
and (c) superconduct as a result of this hybridization. There remain many unresolved 
questions regarding the role of the chains, however. First, there are practical issues: 
neither the strength nor momentum- dependence of the chain-plane coupling have been 
established with any degree of certainty. Second, there are questions of fundamental 
interest: what are the properties of a weakly correlated ID metal in close contact 
with a strongly-correlated electron liquid. Finally, there is the question of how the 
chains manifest themselves in various experiments. In some cases, for example transport 
anisotropy, the role played by the chains is intuitively obvious. However, in many cases, 
the effects of the chains are not a priori obvious, as in the case of the vortex core 
contraction mentioned above. 

In this work, we examine the effects of chains on the local density of states 
(LDOS) p(t,u) in the vortex state of YBa 2 Cus07. Experiments have found that the 
electronic structure of vortices in YBa2Cu30 7 _,5 is different from other HTS[23] and 
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Figure 1. Electronic structure near the vortex core for a single tetragonal layer with 
a d-wave superconducting order parameter, (a) Magnitude of the current density |j(r) 
as a function of position. The figure shows a pair of vortices belonging to a square 
lattice whose axes are rotated 45° relative to the crystalline axes (which are parallel to 
the figure axes), (b) Local density of states at a series of positions along the line drawn 
through the lower vortex core in (a). Note that the nondispersing peak at u> = —0.9 
is a van Hove singularity associated with the band structure. Far from the vortex 
core, the d-wave order parameter is |A| = 0.38. (c) The LDOS p(r,Ep) is shown as a 
function of position near a single vortex core. Note the logarithmic scale. Here and 
throughout the Fermi energy is Bp = 0. 

from theoretical calculations [241125] . Theory predicts that, for a single superconducting 
layer with a c?-wave order parameter, there should be a peak in the LDOS at the vortex 
core centre (a "zero bias conduction peak" or ZBCP), and that as one moves away from 
the core, the peak splits and disperses towards higher energies, eventually merging with 
the coherence peaks far from the vortex core. This behaviour is illustrated in Figure [U 
where the results of a calculation performed for a single-layer version of the model used 
in this work are shown. 

Experimentally, the situation is rather different. The ZBCP at the vortex core 
is absent, and a pair of finite-energy peaks at uj ~ ±5 meV are present instead. 
Furthermore, unlike the calculations, these peaks do not shift significantly as one moves 
away from the vortex core. It is possible that these peaks have the same origin as peaks 
measured in zero-field tunneling experiments, but this is not established. 

In this work we ask two questions: (i) what are the signatures of the chains in the 
vortex structure within the proximity model and (ii) to what extent can the proximity 
model explain the existing experiments? 
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2. Method 

Since the technical details of the calculation have been given in [19], we only give a 
summary of the main ideas here. The plane-chain model is made up of a single bilayer 
consisting of a tetragonal 2D plane layer (representing a CuC>2 layer) and a layer of ID 
chains (representing a CuO layer) . The layers are connected by single electron hopping, 
so that the plane and chain bands hybridize to form new bands. These bands have 
either predominantly plane or chain character, and we refer to them as plane-derived or 
chain-derived respectively. The chain layer is represented by a tight binding Hamiltonian 
with a single orbital per unit cell. In band structure calculations, this orbital is a linear 
combination of the chain Cu^ 2 _ 2 and Py orbitals.pQ The Cu0 2 band structure is also 
represented by a tight binding model with a single orbital per unit cell. In this case, 
the orbital represents a singlet state involving the Cud z2 _ 2 and P(xy) orbitals.[26j. 

There is a pairing interaction in the plane layer leading to (i-wave superconductivity, 
but the chain layer is intrinsically normal. Nonetheless, the chain superconducts because 
of the the proximity effect; single electron hopping between the plane and chain layers 
induces a gap in the chains. One important feature of the proximity model is that the 
induced gap in the chains does not share the ci-wave symmetry of the plane layer [15] . 
This is because the chain dispersion is quasi-lD, and the k-dependence of the excitation 
gap depends on the symmetry of both the order parameter and the underlying band. 

We define coordinates x, y, and z which are aligned with the crystalline axes 
a, b, and c respectively, and consider a magnetic field aligned with the c-axis. The 
screening currents, therefore, circulate within the plane and chain layers. It is for this 
configuration that the vortex core contraction is seen in /zSR experiments. |16j and that 
STM experiments are performed. 

The Hamiltonian is 

H = H 1 + H 2 + H L . (1) 

where Hi is the Hamiltonian for the isolated plane, if 2 the Hamiltonian for the isolated 
chains, and the single-electron hopping term that couples the two layers. For 
comparison, we also consider a single-layer model described by H\ alone. We have 

ijcr ij 

+A*. Cli (r,-)c 1T fo)], (2) 

where ci CT (rj) is the annihilation operator for an electron in the plane on site i with spin a, 
and position = (xi, yi), tuj are hopping matrix elements, and A^- are superconducting 
pair energies. The subscripts "1" and "2" refer to the plane and chain layers respectively. 
The hopping matrix element tuj between sites i and j includes the effects of the magnetic 
field via 

t\ij = tuj exp 



e f Tl 
i— / dr ■ A(r) 
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where tuj are the zero-field matrix elements, A(rj) = — B^yix. is the static magnetic 
vector potential, Bq is the uniform applied magnetic field, and a = eBo/hc. 

We take a square tight-binding lattice with hopping up to second-nearest neigbours. 
The matrix elements in zero magnetic field are tm = ti j0 , = ^i,nn, an d = 

ii.nnn, where and {(i,j)) refer to nearest and next-nearest neighbours respectively. 
When B = 0, the dispersion is e^k) = t lfi + 2t liIin (cos k x + cos k y ) + 4t linnn cos k x cos k y . 

The order parameter for superconductivity is denoted Ajj, where i and j denote 
sites on the lattice. The order parameter is determined self-consistently using a nearest- 
neighbour attractive interaction of magnitude V. Then, 



where t2u = t 2 ,o and t 2i j = t 2jim for i and j nearest-neighbour sites belonging to the 
same chain. The matrix elements are unchanged by the magnetic field because of the 
choice of gauge. When B = 0, the chain dispersion is e 2 (k) = t 2 ,o + 2t 2 ^ m cosk y . 
The term describing the interlayer hopping is 



which mixes the chain and plane wavefunctions. 

The model parameters used in this work are {ti,o, ^i,nn, ^i,nnn, ^2,o^2,nn, t±} = 
{1.0,-1.0,0.45,2.4,-2.0,0.6}, and the pairing interaction is V = 1.3. With these 
definitions, we have taken the magnitude of the nearest-neighbour hopping (i.e. |£i, im |) 
as the scale of energy. This will be our energy scale throughout this work, and a rough 
comparison to experiments may be made by taking ti iim ~ 100 meV. However, we 
emphasize that quantitative comparisons to experiments are not possible because the 
model parameters are chosen for numerical convenience (i.e. such that the energy scale 
i? s associated with chain superconductivity is easily resolved), rather than to reproduce 
the YBa2Cus07_5 band structure accurately. 

There is a quasi-periodicity which allows us to define an L x x L y magnetic supercell 
containing iV = L x L y ja\ atomic lattice sites (a is the lattice constant) and enclosing 
two superconducting flux quanta (i.e. two vortices), where the superconducting flux 
quantum is $0 = hc/2e. The magnetic field is therefore B = 2Q /L x L y . In order 
to obtain low magnetic fields, we need large values of L x and L y . Unless otherwise 
stated, all results shown in this paper are for a 60 x 60 x A^-site supercell (N z = 1,2) 
corresponding to B = 2$ /3600. For this case, we sum over 5 2 = 25 K- vectors, which 




The d-wave component, defined by A(r,) = l) Vi ~ V] ^ij-< is the dominant component 

of the order parameter. 

The Hamiltonian for the chain layer is 

H2 = Y] t 2 i j cl a (r i )c 2 a(r j ) (5) 





(6) 
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is formally equivalent to studying a system of 50 vortices on a 300 x 300 x 2 site lattice. 
Taking a unit cell size of a$ = 4 A, this corresponds to B — 7 T. While this is a strong 
field for YBa 2 Cu307_5, we have chosen the model parameters such that this is in the 
low-field limit of the proximity model (i.e. chain superconductivity is not quenched at 
this field). 

The details of the transformation to Bloch states are given in [T9j. Here, we 
simply mention that our self-consistent calculations proceed in several steps. First, we 
diagonalize the Hamiltonian to generate a set of eigenenergies E^k arid wavefunctions 
^a,K{n, c, r) where K are the supercell wavevectors, a G (1, 27V) (the factor of two is for 
spin) is a quantum number labelling the eigenenergies for each K, n — 1, 2 refers to the 
layer, and r to sites within each layer. Note that a Bogolyubov transformation has been 
made, such that the spin index refers to spin-up electrons with wavevector K (a =f) 
or spin-down holes with wavevector — K (cr =!). The second step is to calculate the 
Fourier-transformed order parameter Ajj(K) from the wavefunctions. This new order 
parameter is inserted into the Hamiltonian, and a new set of eigenfunctions and energies 
is calculated. The iterative process terminates when the largest difference between the 
input and output values of Ajj(K) is less than 10~ 3 . 

We present results for two observables in this work. The LDOS at energy u in layer 

n is 




+\^ a>K {n,i,r)\ 2 S(u + E a>K )] , (7) 

and the total density of states (DOS) is p(u) = N' 1 J2 n S r Pn( r , where N is the 
total number of sites in the lattice. The second observable is the 2D current density in 
layer n, denoted by j„(r) = (j n>x (r), j n , y (r)) with 

. 2N 
U K a=l 

x lYl &.ii( K )*a,K( n . T, ri)* Q)K (n, T,r; ± a x) 
1 ± 

-%(-K)* Q)K (n, |, ri)*; iK (n, |, r, ± a x)] , (8) 

and a similar expression for j yn (ri). Note that the layer index n is omitted when we 
present quantities for the single-layer model. 

3. Results and Discussion 

The LDOS is shown in Figure [1] for a model consisting of a single tetragonal plane. 
This case is well-studied, and we include it here as a point of comparison for later 
calculations. In this figure, one sees a ZBCP at the vortex core. Away from the vortex 
core, the ZBCP splits into two peaks which disperse away from the Fermi energy Ep 
as the distance to the core increases, eventually merging with the coherence peaks at 
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Figure 2. Superconducting density of states, (a) Spatially-averaged density of states 
in zero magnetic field for the bilayer model. The superconducting coherence peaks for 
the 2D plane are split by plane-chain coupling (A), while there are two distinct energy 
scales (B and C) associated with the induced gap in the chains. The inset shows a 
sketch of the chain (red) and plane (blue) Fermi surfaces. The portions of the chain 
Fermi surface contributing to the DOS at the energies B and C are indicated, (b) 
The effect of magnetic field on the spatially-averaged DOS. The field strengths are 
B = 2$ /60 2 (B60) and B = 2$ /30 2 (B30). The latter results are calculated using a 
30 x 30 x 2 lattice with 15 2 supercell K-points. 

uj = ±El. This is the behaviour predicted for isolated vortices in both s-|27j and d- 
wave[25 | 1241 [28] superconductors. For an isolated vortex in an s-wave superconductor, 
quasiparticle eigenstates have a well-defined angular momentum, and it was shown that 
peaks at different distances from the core correspond to different angular momentum 
eigenstates. Sufficiently near the core, both the peak energy and distance to the core are 
linearly proportional to the angular momentum. [29] As we will show, this relationship 
between angular momentum and the peak energy is relevant to the electronic structure 
of CuO chains in the vortex state. 

We remark that, although dispersing peaks have been observed experimentally in 
the s-wave superconductor NbSe 2 |30j, they have not been seen in either Yl^CusOr^ 
or Bi2Sr2CaCu20s423j A number of strong-correlation-related mechanisms have been 
proposed for the absence of the ZBCP[23], and it is known that disorder can suppress 
the ZBCP [23]. Furthermore, it is possible that nontrivial tunneling matrix elements can 
significantly alter the LDOS. [31] At present, the reason for the lack of a ZBCP in the 
HTS has not been resolved. 

The spatial dependence of the LDOS at the Fermi energy, p(r, ep), is plotted in 
Figure (He). A logarithmic plot is used to emphasize the long-range tails extending from 
the vortex core. These are a signature of BCS-like ci-wave superconductivity [25] and, 
like the ZBCP, have not been seen experimentally. Again, the reason for this is not 
clear. 
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In Figure [2J we compare the spatially-averaged density of states for three different 
cases: the bilayer model in zero magnetic field, and the single-layer and bilayer 
models for B ^ 0. This figure shows that the DOS for the bilayer model has a 
relatively complicated structure when B = 0, with four distinct energy scales (see 
[T5] for a complete discussion). First, the coherence peaks associated with plane 
superconductivity are split by chain-plane coupling (A). A similar splitting has been 
observed in STM experiments [22J. We identify the lower-energy of the two coherence 
peaks as the "large gap" E^. Because the chains act as pair-breakers, £7l is slightly 
smaller than the gap in the single-layer model. Next, there are two gap energy scales 
associated with induced superconductivity in the chains. The larger of these (B) comes 
from the region of the chain-derived Fermi surface nearest (jr, 0) in the Brillouin zone. 
Here, the energetic proximity of the chain and plane bands leads to a large hybridization 
of their wavefunctions and a relatively large induced gap. We believe that this energy 
scale, denoted E s , is the small gap that has been seen in STM experiments in zero field. 
The smallest energy scale (C) is associated with superconductivity on the portion of the 
chain Fermi surface nearest the Brillouin zone centre. This gap is wiped out by small 
amounts of disorder or thermal broadening, and the proximity model therefore predicts 
a finite residual density of states at the Fermi energy, in accordance with tunneling 
experiments. [201I2TII22] 

The finite- field DOS is shown in Figure [2(b). The field washes out the distinct 
energy scales; however, E s shows up in the field-dependence of the DOS. Compared with 
the single-layer model, p{ef) for the bilayer model is a strong function of B for B < B*, 
where B* ~ E*\£5\ is the crossover field mentioned above. For B > B*, p(ep) varies 
at a similar rate with B in both the single-layer and bilayer models. Similar behaviour 
has been predicted in an anisotropic single-band model for YBa9CuqQ7_A[32]. In the 
bilayer model studied here, B* is the magnetic field above which the vortex cores in 
the chains begin to overlap. From the perspective of the chains, there are therefore 
distinct low- {B < B*) and high- (B > B*) field regimes. Experimentally, B* « 1 T in 
YBa 2 Cu 3 06.95- Because we have chosen a large value for the chain-plane coupling, the 
data shown here for the 60 x 60 x 2 lattices belong to the regime B < B* . 

The LDOS for the bilayer is shown along a series of cuts through the vortex along 
the (1, 0) direction in Figure El the (1, 1) direction in Figure EJ and the (0, 1) direction in 
Figure [5j The cuts are illustrated in Figure [3]^d). There are three sources of anisotropy 
which contribute to differences in p n (r,uj) along the three directions. First, the d- 
wave order parameter vanishes for quasiparticles travelling in the (1, ±1) directions and 
obtains its maximum along the (1, 0) and (0, 1) directions. Second, the shortest path 
between vortices lies along the (1, ±1) directions. Third, the chains run parallel to the 
(0, 1) direction. We can separate the third factor from the first two by comparing our 
bilayer results with those for the single-layer model. 

First, we consider the cut along the (1,0) direction (Figure [3]). The LDOS in the 
plane layer exhibits one notable difference from the single layer case shown in Figure [U 
the ZBCP splits into a quartet, instead of a pair, of dispersing peaks as one moves away 
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Figure 3. Local density of states in the bilayer model, (a) The LDOS pi(r, uS) in the 
2D plane at positions between the vortex core (top curve) and the edge of the vortex 
(bottom curve). Curves are offset for clarity. Positions are taken along the (1,0) cut 
shown in (d). (b) p2(r,u>) in the ID chain layer at identical positions to (a), (c) Two 
of the LDOS curves from (a) are compared, showing the dispersion of the inner (I) and 
outer (O) peaks. Curve colours are the same as in (a), (d) Spatial dependence of the 
current density |ji(r)| in the plane layer for a pair of vortices. The lines indicate the 
cuts along which LDOS curves are displayed in this figure [(1, 0)], Figured] [(1, 1)] and 
Figure [5] [(0, 1)]. Note that the (0, 1) and (1,0) cuts are parallel and perpendicular to 
the chains respectively. 

from the vortex core. The outer peaks asymptotically approach E^ at large distances 
from the vortex core, while the inner peaks asymptotically approach E s . It appears as if 
the presence of two superconducting energy scales leads to the formation of two distinct 
sets of dispersing quasi-bound resonances. 

The LDOS curves shown in Figure [3](b) come from different chains at increasing 
distances from the vortex core. As in the plane, there is a ZBCP at the vortex core. The 
splitting of the ZBCP that is readily apparent in the plane layer is also present, although 
difficult to resolve, in the chain layer. The LDOS curves far from the vortex core closely 
resemble the DOS for the chains when B = 0: there is a pronounced suppression of the 
DOS for \u\ < E s , a large residual DOS at u = ep, and there are weak remnants of the 
coherence peaks at E^ coming from the superconductivity in the planes. 

The LDOS is shown along the (1, 1) direction in Figure HI In both the single-layer 
and bilayer models, the peaks disperse more quickly along the (1, 1) direction than along 
the (1, 0) direction. Furthermore, in the case of the bilayer model, the inner peaks are 
weak and difficult to resolve. 

When the cut is taken along the (0, 1) direction (Figure [5]), the single layer model 
gives the same LDOS curves as in the (1,0) direction, as required by the fourfold 
symmetry of the lattice. The bilayer model, however, is qualitatively different along 
the (0, 1) and (1, 0) cuts. In the plane layer, there are a pair of dispersing peaks 
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Figure 4. Local density of states in the bilayer model for the (a) plane and (b) chain 
layers. Cuts are taken along the (1,1) direction, illustrated in Figure [D^d) , with the 
vortex core corresponding to the top curve. The LDOS is also shown in (c) for the 
single-layer model. 
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Figure 5. Local density of states in the bilayer model for the (a) plane and (b) chain 
layers. Cuts are taken along the (0,1) direction, illustrated in Figure WiA) , with the 
vortex core corresponding to the top curve. The LDOS is also shown in (c) for the 
single-layer model. 
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Figure 6. LDOS maps for two vortices on a 60 x 60 x N z lattice. Figures show 
p n (r,u>) as a function of r for the single-layer (N z — 1; top row), and bilayer model 
(N z = 2; middle and bottom rows). Each column corresponds to a different value of 
ui, indicated at the top of the column. 

that asymptotically approach ±Ei at large distances from the vortex core. These 
peaks correspond to the outer peaks found along the (1, 0) direction. In contrast, the 
inner peaks that are seen along the (1, 0) direction are absent here. Instead, there 
is a pronounced ZBCP that extends ~ 10 unit cells from the vortex core along the 
(0, 1) direction. The ZBCP is particularly pronounced in the chain layer, where extends 
a longer distance than in the plane layer. Closer examination reveals a very slight 
splitting of the peak at large distances. For reasons discussed below, we believe that 
this splitting occurs because chain-plane coupling makes the chain-derived states weakly 
2D. 

For comparison, the spatial maps of p n (r,uj) are shown at three different values of 
u in Figure El This figure provides a different perspective on the results presented in 
the preceding discussion. We begin with the plots for uj = Sp. The single-layer model 
exhibits a fourfold symmetry, as expected. The ZBCP appears as the bright cross in the 
middle of each vortex core, and disappears away from core, where the peaks disperse to 
higher energies. In the bilayer model, we can clearly see that the LDOS in the plane, 
Pi(r, £p), is orthorhomically distorted and has long tails along the chain direction. In 
the chain layer, the vortex core is even more anisotropic; the extent of the core states 
along the a-axis is similar to that in the plane layer, but the extent along the 6-axis is 
determined by the coherence length in the chains, £ c = HvF tC /7rE s m 10a , where vp tC is 
the chain Fermi velocity. 

At the two larger values of u, the LDOS shows a vortex core that appears to expand 
with increasing u. This is because the apparent border of the core is determined by the 
positions of the dispersing peaks at that value of uj. As discussed earlier, the position and 
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energy of these peaks are approximately linearly related. The most noticeable difference 
between the single-layer and bilayer models is the orthorhombic distortion of the vortex 
core in the bilayer model, not present in the single-layer calculation. The single-layer 
model exhibits the fourfold "star" shape, first described in [25], while the LDOS in the 
bilayer plane is stretched along the a-axis, perpendicular to the chains. We note that 
STM experiments in YBa2Cu30y_5 found a similar stretching of the cores along the 
a-axis, with a ratio of ~ 1.5 between major and minor axes. [23] Although our results 
are suggestive, it is premature to declare the chains to be the underlying mechanism for 
this anisotropy, especially given the number of unresolved discrepancies between theory 
and experiment. We return to this point below. 

We can construct a qualitative argument for the anisotropy between the (1, 0) 
and (0, 1) directions in the dispersion of the inner peaks. The fact that these peaks 
asymptotically approach E s at large distances from the vortex core suggests that they 
are closely connected to superconductivity in the chain layer. We therefore consider the 
behaviour of quasiparticles travelling at the chain Fermi velocity vy >c within a chain a 
perpendicular distance r±_ from the vortex core. Since the direction along which the 
quasiparticle travels does not change as it moves along the chain, the quasiparticle 
angular momentum is conserved and is L c = rxvp >c . This is obviously an approximate 
statement since a quasiparticle, as a result of chain-plane hybridization, spends a fraction 
of its time in the plane layer where its trajectory is altered by the magnetic field. 
The linear relationship between r± and L c is therefore also approximate. However, it 
is sufficient to explain the anisotropy between (1, 0) and (0, 1) directions. We recall 
that in conventional superconductors there is a linear relationship between the angular 
momentum of a quasiparticle eigenstate and its energy. [29l 127] For the (1, 0) cut, r± is 
equal to the distance r from the vortex core, L c grows linearly with r, and the inner peaks 
are expected (by this argument) to disperse linearly with r. This prediction is consistent 
with the numerical results shown in Figure [3j On the other hand, r± = everywhere 
along the (0, 1) cut, since that cut follows along the chain that passes through the vortex 
core. We therefore have L c = everywhere along this cut, and the inner peaks remain 
at uj = Ep. Since the single-electron hopping between the chain and plane layers makes 
the chain-derived states weakly 2D, L c is only approximately conserved along the length 
of the chain. We therefore expect a very weak dispersion of the peak energies along the 
(0, 1) direction. Again, this simple argument is consistent with numerical results, shown 
in Figure El 

We finish this section with a brief discussion of the relevance of this work to 
experiments. There is a large body of literature attempting to explain STM experiments 
in HTS, both in zero magnetic field and in finite fields [23]. Existing theories are 
successful at reproducing broad qualitative features of the DOS spectrum, but generally 
fail to correctly predict quantitative details such as the spatial dependence of the 
LDOS near isolated impurities or vortex cores. The same appears to be true of the 
work discussed here: the proximity model is consistent with experiments indicating 
the existence of a small superconducting gap scale E s « 5 meV but, as with previous 
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calculations, there is a discrepancy between the experiments and theory regarding the 
spatial dependence of the LDOS. Whether the failure to explain experiments is an 
experimental issue, or is due to unknown tunneling matrix element effects, surface states, 
or a failure of the basic mean-field model, is unclear at this stage. 

Because of these uncertainties, we seek simple qualitative predictions stemming 
from the proximity model that might reveal something about the electronic structure 
of the vortex cores. We believe that one qualitative feature predicted by the proximity 
model, the existence of an anistropy in the low energy LDOS between the (1,0) and 
(0, 1) directions, should hold regardless of the details of the model. An experiment to 
measure such an anisotropy would establish whether E s does, in fact, originate with the 
chains. In order to address this question, it will be necessary to perform experiments at 
strong fields where the vortex lattice is square. 

4. Conclusions 

We have studied the local density of states near a vortex core within a proximity model 
for YBa 2 Cu 3 7 _ < 5. This model incorporates both Cu0 2 planes and CuO chains. The 
proximity model predicts that the CuO chains introduce a set of dispersing peaks in 
the LDOS, in addition to the peaks predicted for a single-layer superconductor. These 
peaks are associated with the small induced gap in the chain layer, and consequently 
have a strongly anisotropic dispersion. At uo = ef, the LDOS is extended along the 
chain direction, but at uj > e f , the vortex cores appear stretched along the a-axis. 
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